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We study a lattice model where particles carrying different masses diffuse, coalesce upon contact, 
and also unit masses adsorb to a site with rate q or desorb from a site with nonzero mass with 
rate p. In the limit p = (without desorption), our model reduces to the well studied Takayasu 
model where the steady-state single site mass distribution has a power law tail P(m) ~ mT T for 
large mass. We show that varying the desorption rate p induces a nonequilibrium phase transition 
in all dimensions. For fixed q, there is a critical p c (q) such that if p < p c (q), the steady state mass 
distribution, P(m) ~ m" T for large m as in the Takayasu case. For p = p c (q), we find P(m) ~ m~ Tc 
where r c is a new exponent, while for p > p c (q), P(m) ~ exp(— m/m*) for large m. The model is 
studied analytically within a mean field theory and numerically in one dimension. 
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I. INTRODUCTION 



Many systems in nature, ranging from reaction-diffusion systems to fluctuating interfaces, exhibit nonequilibrium 
steady states with a wide variety of phases. Of particular interest are the self-organized critical systems where different 
physical quantities have power law distributions in the steady state over a wide region of the parameter space E). 
Self-organized criticality has been studied in a variety of model systems ranging from sandpiles to earthquakes. A 
particularly simple lattice model due to Takayasu, where masses diffuse, aggregate upon contact and adsorb unit 
masses from outside at a constant rate, was shown to exhibit self-organized criticality jl: the steady-state mass 
distribution has a nontrivial power law decay for large mass in all dimensions . This model initially generated a lot 
of attention as it was a simple exactly solvable model of self organized criticality with close connections || to other 
solvable models such as the Scheidegger river model Q, the voter model j^] and the directed abelian sandpile model 
H . Recently there has been a renewed interest in this model as simple variants of the Takayasu model have been 
found useful in modelling the dynamics of a variety of systems including force fluctuations in granular systems such 
as bead packs 0, river networks ||, voting systems [P,|lO|], wealth distributions pL inelastic collisions in granular 
gases Jl3| , the generalized Hammersley process [fL2fl , particle systems in one dimension [ fl4| and various generalized 
mass transport models p5| . 

In the Takayasu model, each site of a lattice has a nonnegative mass variable. Starting from an initial random 
distribution of masses, each mass hops to a nearest neighbour site (chosen at random) and aggregates with the mass 
there with rate 1. In addition, a unit mass is adsorbed at every site with rate q. While the first move tends to 
create big masses via diffusion and aggregation, the second move replenishes the lower end of the mass spectrum. At 
large time t, the mass distribution at any site has the scaling behaviour, P(m,t) ~ mT T f(m/t s ) with S = 1/(2 — r) 
jl6],[l7|]. The interesting point is that even though the average mass per site increases linearly with time, (m) <~ t, 
the mass distribution P(m,t) approaches a time-independent power law distribution P(m) ~ mT T for t — > oo (since 
/(0) ~ O(l)) for any nonzero adsorption rate q. The exponent r is independent of q and is known exactly r = 4/3 
in one dimension and r = 3/2 within mean field theory. 

The steady-state mass distribution in the Takayasu model has the same power law decay for any nonzero adsorption 
rate q and does not undergo any phase transition. In this paper we show that if we introduce an additional process 
of desorption of unit masses with rate p in the Takayasu model (we call this the In-out model), a rich steady-state 
phase diagram emerges in the p — q plane. In particular we show that the system undergoes a nonequilibrium 
phase transition across a phase boundary p c {q)- Nonequilibrium phase transitions between steady states have been 
studied extensively in recent years in a variety of systems. Examples include, amongst others, active-absorbing phase 
transitions in reaction diffusion systems [ fl8f , roughening transitions in fluctuating interfaces [fl9f , phase transitions in 
driven diffusive lattice gas models f20(| , wetting transitions in solid-on-solid models |21|, boundary driven transitions 
in one dimensional asymmetric exclusion processes and Bose-Einstein like condensation in models of aggregation 
and fragmentation |2^,|24|] . However we show below that the mechanism of the phase transition and the associated 
critical properties in the In-out model are very different from those of other models mentioned above. 
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There are quite a few physical systems where our In-out model may find applications. In nature there exist a 
variety of systems ranging from colloids [p5| to polymer gels [^6| where the basic constituents of the system diffuse 
and coalesce upon contact. For example, in a polymer gel the basic constituents are polymers of different sizes which 
diffuse in a solution and when an m-mer comes in contact with an n-mer, they aggregate to form an (m + n)-mer 
p6[ . Similarly during the growth of a thin film on an amorphous substrate (such as Bismuth on Carbon) , clusters or 
islands of atoms can diffuse as a whole and when two of them come closer they coalesce |27) . A zeroth order approach 
to model the dynamics of these systems would be to replace each cluster by a point particle (ignoring its shape) 
carrying a positive mass which indicates its size or number of atoms. When two particles coalesce their masses add 
up. In addition many of these systems are open in the sense that they can exchange basic units with the adjoining 
environment. For example, during the growth of a film on a substrate, single adatoms may adsorb on the substrate 
from the outside vapour or desorb into the vapour from the substrate. We attempt to incorporate these processes on 
a lattice in the In-out model and show that even this simple model has a very rich steady-state phase diagram. We 
had introduced this model in an earlier publication fjl| and some results were briefly mentioned. In this paper we 
present a detailed analysis of the model. 

The paper is organized as follows. In section II, we define the In-out model precisely and summarize the different 
phases and the transitions between them. In section III, we solve the model analytically within mean field theory. 
In section IV, we present the numerical results in one dimension and discuss a scaling theory which provides scaling 
relations between different critical exponents. We conclude in section V with a summary and a discussion of open 
questions. 



II. THE IN-OUT MODEL 



For simplicity we define the In-out model on a one-dimensional lattice with periodic boundary conditions; gener- 
alizations to higher dimensions are straightforward. Each site of a lattice has a nonnegative mass variable > 0. 
Initially each rrii is chosen independently from any well defined distribution. The dynamics proceeds as follows. A 
site i is chosen at random and then one of the following events can occur: 

1. Adsorption: With probability q/(p + q + 1), a single particle is adsorbed at site i; thus rrii — » m t + 1. 

2. Desorption With probability p/(p + q + 1), a single particle detaches from and leaves site i; thus m., — > m., — 1 
provided m., > 1. 

3. Diffusion and Aggregation: With probability l/(p + q + 1), the mass wii at site i moves to a nearest neighbour 
site [either (i — 1) or (i + 1)] chosen at random. If it moves to a site which already has some particles, then the 
total mass just adds up; thus rrii — > and rrii±i — » mm + rrii. 

If the site chosen is empty, only adsorption can occur with probability q/(p + q + 1). 

The In-out model has only two parameters p and q. The question we would like to address is: For given p and q, 
what is the single site mass distribution P(m) in the steady state? Note that in the limit p = (i.e., without the 
desorption process) our model reduces to the Takayasu model mentioned in the introduction. 

While the Takayasu model (zero desorption, p = 0) does not have a phase transition in the steady state, we find 
that introducing a nonzero desorption rate p induces a rich steady state behaviour in the p — q plane. In fact we find 
that there is a critical line p c {q) in the p — q plane. For fixed q, if we increase p from 0, we find that for all p < p c (q), 
the steady state mass distribution has the same large m behaviour as in the Takayasu case, i.e., P(m) ~ m~ T 
where the exponent r is the Takayasu exponent and is independent of q. Thus the Takayasu phase is stable upto 
p c . For p — p c (q), we find the steady state mass distribution still decays algebraically for large to, P{m) ~ mT Ta 
but with a new critical exponent r c which is bigger than the Takayasu exponent r. For p > p c (q), we find that 
P(m) ~ exp(— to/to*) for large m where to* is a characteristic mass that diverges if one approaches p c {q) from the 
p > p c {q) side. The critical exponent t c is the same everywhere on the critical line p c {q)- This phase transition occurs 
in all spatial dimensions including d = 1. 

It is easy to write down an exact evolution equation for the mean mass (m)(t) per site. Since the diffusion and 
aggregation move does not change the total mass, the only contributions to the time evolution of (m) come from the 
adsorption and desorption processes. It is then evident that, 

——=q-ps{t) (1) 
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where s(t) is the probability that a site is occupied by a nonzero mass. The first term on the right hand side of the 
above equation clearly indicates the increase in mass per site due to the adsorption of unit mass. The second term 
quantifies the loss in mass per site due to the desorption of unit mass taking into acount the fact that the desorption 
can take place from a site only if the site is occupied by a nonzero mass. Let us fix p and vary q. As long as q < q c {p), 
it turns out that in the long time limit t — ► oo, the two terms on the right hand side of the above equation cancel 
each other and the occupation density reaches the asymptotic time independent value, s = q/p. This indicates that 
the average mass per site, (m) becomes a constant in the long time limit. In fact, we show below that in this phase, 
the steady state mass distribution P{m) ~ exp(— m/m*) for large m with a finite first moment (m). We call this 
phase the "Exponential " phase. However if q > q c (p), the occupation density reaches a steady state value s such 
that s < q/p. As a result in the long time limit, the second term on the right hand side of Eq. (1) fails to cancel 
the first term and the mean mass per site (m)(t) increases linearly with time, (m) ~ (q — ps)t. However, as we show 
below, even though the mean mass diverges in this phase as t — > oo, the mass distribution reaches a steady state, 
P(m) ~ mT T for large m where r is the Takayasu exponent (which is always less than 2 so that the mean mass 
diverges). Hence we call this entire phase the "Takayasu" phase. 



III. MEAN FIELD THEORY 



We first analyze the model exactly within the mean field approximation, ignoring correlations in the occupancy of 
adjacent sites. In that case we can directly write down equations for P(m, t), the probability that any site has a mass 
m at time t. 



dP(m, t) 
dt 

dP(0,t) 

dl 



(1 + p + q + s)P(m,t) +pP(m + l,t) 

qP(m- l,t) + P * P; m > 1 (2) 
(q + s)P(0,t)+ P P(l,t) + s(t). (3) 



Here P * P = Ylm'=x P( m ' ,t)P{ m — m',t) is a convolution term that describes the coalescence of two masses and 
s(t) = Yl m —x P(m,t) denotes the probability that a site is occupied by a nonzero mass. 

The above equations enumerate the possible ways in which the mass at a site might change. The first term in Eq. 
(|J) is the "loss" term that accounts for the probability that a mass m might move as a whole or desorb or adsorb 
a unit mass, or a mass from the neighbouring site might move to the site in consideration. In this last case, the 
probability of occupation of the neighbouring site, s(t) multiplies P{m, t) within the mean-field approximation where 
one neglects the spatial correlations in the occupation probabilities of neighbouring sites. The remaining three terms 
in Eq. (^|) are the "gain" terms enumerating the number of ways that a site with mass m' ^ m can gain or lose mass 
to make the final mass m. The second equation Eq. (|3|) is a similar enumeration of the possibilities for loss and gain 
of empty sites. 

To solve the equations, we compute the generating function, Q(z, t) = J2m=i P{ m -> i) 2 ™ horn Eq. (|2|) and set 
dQ/dt — in the steady state. We also need to use Eq. (|^) to write P(l, t) in terms of s(t). This gives us a quadratic 
equation for Q in the steady state. Choosing the root that corresponds to Q(z = 0) = 0, we find 



2zQ(z) = p(z - 1) + qz(l - z) + 2sz - y/(z - l)A(z). (4) 

where 

AO) = p 2 (z - 1) + q 2 z 2 (z - 1) - 2pqz(z - 1) 

-4qz(z- sp/q). (5) 

Note that the occupation density s in the above expression of Q(z) is yet to be determined. The steady state mass 
distribution P(m) can be formally obtained from Q(z) in Eq. (Q) by evaluating the Cauchy integral, 



v ; 2m J c z m+1 



dz (6) 



over a contour C Q encircling the origin in the complex plane. This expression for P(m) however will contain the yet 
to be determined unknown quantity s. In fact, determining s is the most nontrivial part of the mean field calculation 
as we show below. 
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In order to extract the large-m behaviour of P{m) from Eq. (g), one needs to deform the contour C Q so that it goes 
around the branch cut singularities of the function Q[z), From Eq. (H), it is evident that such singularities occur at 
z = 1 and also at the roots of A(z) = where A(z) is given by Eq. ([|). Since A(z) is a cubic polynomial in z, it has 
three roots Zi, z<i and Z3, each of which can be determined in terms of the unknown quantity s. 

We now analyse the large m behaviour of P(m) in different regions of the p — q plane. Let us fix the value of p and 
increase q from 0. A similar analysis can be carried out for fixed q as a function of p. As we increase q from 0, we 
encounter the following three regimes, 

(i) For small q (with a fixed p), we first assume that the mean mass (to) reaches a time- independent constant as 
t — > 00. This assumption will be justified a posteriori. Then from Eq. (1), it follows that the occupation density 
also reaches a steady state value, s = q/p. Substituting this in the expression for A (2) in Eq. (||), the three roots 
of A(z) = are z\ = 1 and 2:2,3 = {p + 2 =F + T)/q. Then from Eq. (Q), it follows that the only branch cut 
singularities of Q{z) are at zi and Z3 with Z3 > zi > 1 for small q. Therefore the branch cut at z^ essentially controls 
the large-m behaviour of P(m) when the contour in Eq. @ is deformed and by analysing the integral around this 
cut we find that for large m, 

P{m) ~ exp(—m/m*)/m 3 / 2 (7) 

with to* = \J\nz%. Since Pirn) decays exponentially in this phase, (to) is also finite and nonzero thus justifying the 
assumption made in the begining. In this phase the unknown function s is therefore exactly given as s = q/p. Note 
however that this analysis is valid as long as Z2 > 1 and the characteristic mass to* diverges as z-i approaches 1 from 
above. 

(ii) As the value of q is increased (for fixed p) the roots Z2 and Z3 decrease, until at a critical value q c (j>), the value 
of Z2 just reaches unity. The double root [zi and z%) at z = 1 of A(z) then leads to a branch cut singularity of order 
3/2 in Q(z) in Eq. (||), which in turn implies 

P(to)~to~ 5/2 . (8) 

This power law decay characterizes the critical point and the condition 22 = 1 determines the locus of the critical line 
in the p — q plane, 

?c(p) =p+2-2y/p+l. (9) 

The value of s is given exactly by s = q c /p. 

(iii) As q is increased further (q > q c {p)) for fixed p, the mean mass per site (to) does not reach a time-independent 
value in the steady state, but increases indefinitely with time. Consequently we cannot use the relation s = q/p 
anymore. However, P(m) reaches a time-independent distribution. So the question is what is the selection principle 
that determines the unknown function s in this regime? 

Note that at q — q c {p), the two roots Z\ and Z2 of A(z) = coincided, Z\ = z% = 1 and Z3 > 1. As q increases 
further, since we do not know what s is a priori, the exact locations of the three roots of A(z) = in the complex 
plane are also unknown. However since A(z) is a polynomial with real coefficients, if z is a root of A(z) = 0, so must 
be its complex conjugate z*. Thus as q increases beyond q c (p), there are two possibilities. The first possibility is that 
all the three roots of A(z) = are real and distinct. But in that case, as q increases slightly beyond q c (p), at least 
one of them must become less than 1. This however would lead to an exponential growth of P(m) for large m and 
hence is ruled out. The second and only possibility is that one of the three roots must be real while the other two are 
complex conjugates of each other, i.e., A(z) = (z — z c )[z — z c *){z — Z3) where Z3 is real and z c in general is complex 
with its real part less than 1. However, if the imaginary part of z c is nonzero, this again can be shown to lead to an 
exponential divergence of P(m) for large m. Therefore, we are led to the conclusion that z c must be real and thus 
A(z) = must have double roots at z c , i.e., A(z) = (z — z c ) 2 (z — Z3) with z c real. In summary we conclude that 
for q > q c (p), Z3 remains > 1 and the two roots z± = zi = z c continues to be coincident and real but the common 
value z c decreases below 1 as q increases beyond q > q c (p)- This nontrivial 'root sticking' condition determines the 
unknown quantity s for q > q c (p). This condition of double roots can be easily implemented by demanding the two 
conditions, A(z c ) = and A'(z c ) = where A' = dA(z)/dz. Also using the relation A(z) = (z — z c ) 2 (z — Z3) in Eq. 

, we find that the lowest branch cut singularity of Q(z) is at z = 1. This order 1/2 singularity then leads to the 
following asymptotic behaviour of P(m), 
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P(m) ~ toT 3/2 . 



(10) 



Thus this entire phase, q > q c {p) is characterized by the same power-law decay of P(m) as in the mean field Takayasu 
model which, as mentioned earlier, corresponds to the zero-desorption (p = 0) limit of our model. 

As mentioned above, the 'root-sticking' condition also determines quite non-trivially the occupation density s for 
q > q c {p) for fixed p. To determine s explicitly for q > q c {p) using this condition, let us fix p = 1 for simplicity even 
though the calculation can be carried out for any arbitrary p. From Eq. (|J), we find q c = 3 — 2 y/2 for p = 1. We first 
substitute the expression for A(z) from Eq. (||) in the 'root-sticking' conditions, A(z c ) = and A'(z c ) — 0. We then 
eliminate z c from these equations and find s(q) for q > 3 — 2>/2 as the only positive root of the cubic equation 

16s 3 - (q 2 - 12q + 24)s 2 - {q 3 + 5q 2 + 57q + 15) s 

+ (q 3 + 5q 2 + 39q-2) = 0. (11) 

Thus we can determine the unknown quantity s exactly everywhere in the p — q plane. In Fig. 1, we plot the function 
s(q) for fixed p = 1. For q < q c = 3 — 2y2 we have s(q) — q and for q > q c = 3 — 2v2, s(q) is given by the real 
positive root of the cubic equation in Eq. (jlT|). 

Note that for fixed p, if q > q c {p), the steady state value of s(q) (as determined from the 'root sticking' conditions) 
is less than q/p and hence from Eq. (1), we find that the mean mass per site increases linearly with time, (m) ~ vt 
for large t. If one interprets the mass profile as the height of an interface (see Section V) then for q < q c (p), the 
average "height" of the interface becomes a constant as t — ► oo, while for q > q c (p), the average "height" (m) increases 
linearly with velocity v. The 'velocity' v defined more precisely as v — limt^oo ^j^- is for q < q c (p) and nonzero 
for q > q c {p)- For q slightly bigger than q c {p), v ~ [q — q c {p)] v where y is a critical exponent independent of p. For 
example for p = 1, we find from Eq. (pT|), v « (q — q c (l)) 2 / (6\/2 — 8) indicating that y = 2 within mean field theory. 



IV. NUMERICAL RESULTS IN ONE DIMENSION AND SCALING THEORY 

Having completed the mean field calculations we now turn to one dimension. While the Takayasu model (p — 0) is 
exactly solvable in d = 1 Q], the same technique unfortunately does not work for p > 0. Hence for nonzero p, we had 
to resort to numerical simulations in d = 1. The qualitative predictions of mean field theory namely the existence 
of a power-law (Takayasu) phase (P(m) ~ mT TT ) and a phase with exponential mass distribution, with a different 
critical behaviour at the transition (P(m) ~ m~ Tc ), are found to hold in 1-d as well. Figure 2 shows the results of 
numerical simulations for the phase diagram along with the mean-field prediction (Eq. |^) and Figure 3 displays the 
numerical data for the decay of the mass distribution P(m) in the two phases and at the transition point. The values 
obtained, r = 4/3 (same as the exactly solvable p = case) and r c ~ 1.833, are quite different from their mean-field 
values t — 3/2 and r c = 5/2, reflecting the effects of correlations between masses at different sites. 

If the phase boundary is crossed by increasing q for fixed p, the Takayasu phase is obtained for q > q c . As a function 
of the small deviation q = q — q c and large time t, the mass distribution P(m, q, t) is expected to display a scaling 
form for large to, 

P( m> g >t )„_L y(mg*£J) (12) 

in terms of three unknown exponents 0, a, r c and the two variable scaling function Y. All other exponents then can 
be related to these three exponents via scaling relations. We give some examples below. 

(a) Consider q > and t — > oo limit. Then P(m,q) ~ ■^^Y{mq^' 1 Q). But we know that for q > 0, in the steady 
state, P(m,q) ~ mT T where r is the known Takayasu exponent. This forces the scaling function Y(x,0) ~ x 7 for 
large x such that P(m,q) ~ g^ 7 / l m Ta ~ 1 ', indicating 7 = r c — r. 

(b) Consider again q > and finite but large t. The mean mass per site, (m) = J mP(m, q, t)dm ~ (pt where y is 
the 'velocity' exponent. Using the scaling form of P, we find, y = 4>[1 — a(2 — r c )]/a. 

(c) Next we consider the critical point, q = 0. Using the scaling form, we find that the mean mass, (to) ~ p for 
large t where £ = a(2 — r c ) provided r c < 2. If r c > 2 (as in mean field theory), £ = 0. Also, the root mean square 
mass fluctuations at the critical point, a = y ((™ ( m )) 2 ) ~ yM ~ ^ f° r l ar g e * with /3 = a(3 — r c )/2. Note 
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that for large t, (m 2 ) >> (m) indicating that fluctuations grow faster than the mean as time increases. 

Within mean field theory, by analysing P(m) explicitly for q > 0, we find P(m,q) ~ qjvr?! 2 and also r c = 5/2. 
From (a) above, this immediately gives, 70 = 1 and 7 = 1 indicating ef> = 1. Also, we had shown before that the 
velocity exponent y — 2 exactly within mean field theory. Using y — 2, t c = 5/2 and (f> = 1 in (b) of the previous 
paragraph, we get a = 2/3. Since r c = 5/2 > 2, we note from (c) that £ = 0. Also we find the fluctuation exponent 
(3=1/6 from the scaling relation in (c). Thus within mean field theory, we find 

i vn 

P(m,q,t) ~ —^Y(mq,-^). (13) 

We have determined the corresponding exponents in d — 1 numerically. The critical exponent t c ~ 1.83 has already 
been mentioned (see Fig 3). In Fig 4, we plot the velocity v as a function of q for fixed p = 2.35. The velocity is zero 
for q < q c ~ 1.0 and increases as a power law, u ~ (q — q c ) v for small q = (q — q c ). We find y ~ 1.47. Note that since 
q c is not known exactly, this exponent is difficult to determine numerically and is subject to large error bars. We also 
find that at the critical point q c ~ 1, the mean mass grows as, (m) ~ ft with £ ~ 0.12. Note the difference from the 
mean field theory where (m) does not grow with time at the critical point (£ = 0). To measure the fluctuations at 
the critical point, we performed finite-size studies of the time-dependent 'width' W 2 (t 7 L) = J2i=i( m i~ < m >) 2 /-^ 
at the critical point, where L is the system size. This is expected to obey the scaling form W ~ ftZ(t/L z ); the value 
of z is expected to be 2 as the movement of masses is diffusive. Figure 5 shows the scaling plot of W/ft versus t/L z 
for four different system sizes L = 16, 32, 64 and 128 at the critical point p w 2.35 for fixed q = 1. We fix z = 2 
and find the best collapse of data for /3 ~ 0.358. These exponent estimates in d — 1 are consistent with the scaling 
relations mentioned in (a)-(c). 



V. SUMMARY AND DISCUSSION 



In this paper we have studied a simple lattice model where masses diffuse and aggregate with rate 1, unit masses 
adsorb at any lattice site with rate q and unit masses desorb from a site (provided the site is occupied by a mass) 
with rate p. For p = (without the desorption process), our model reduces to the well studied Takayasu model where 
the steady-state single site mass distribution has a power law decay, P(m) ~ mT T for large m for any nonzero q. 
We show that varying the desorption rate p induces a nonequilibrium phase transition at a critical value p = p c (q)- 
For p < p c (q)i P(m) ~ mT T for large m as in the Takayasu (p = 0) case. For p = p c (q), P(m) ~ m~ Tc where t c is 
a new exponent and P(m) ~ exp(— m/m*) for p > p c {q)- We have solved the model analytically within the mean 
field theory and calculated all the mean field exponents exactly. In one dimension, we have computed the exponents 
numerically. We have also presented a general scaling theory. 

Our model generalizes the Takayasu model and exhibits a nontrivial phase transition. There was an earlier general- 
ization pjfl of the model where instead of carrying positive masses, the diffusing particles carried charges Q of either 
sign while a random charge /, drawn from an arbitrary distribution, was added with rate q to a lattice site. In this 
"charge" model, the steady-state single site charge distribution P(Q) was found (28) to have a power law tail (as in 
the mass case), P(Q) ~ Q~ T for large positive Q when the mean charge injected was positive, (I) > whereas for 
(7) = 0, P(Q) ~ Q~ Tl for large positive Q. It was shown that the exponent t\ = 5/3 in d — 1 and t\ = 2 within 
mean field theory [M . Though this change of exponent at a critical value (I) = is similar to that in our model 
qualitatively, the exponent r c of the In-out model is very different from that of the "charge" model. This difference 
can be traced back to the mass positivity constraint in the In-out model, i.e., the desorption of an unit mass can take 
place from a lattice site only if the site has a nonzero mass. 

In the In-out model, the total mass is not conserved due to the moves involving adsorption and desorption of 
unit mass. It is interesting to ask what would happen if the desorption of a unit mass from a site were followed by 
adsorption at a neighbouring site, so that the total mass would be conserved in every move. This was investigated 
using a lattice model pj| and earlier, within a rate equation approach [M. In this conserved-mass model too there is 
a phase transition, but of a different character. It was found that there is an exponential phase (at high desorption- 
adsorption rate), separated by a critical line from a phase with a power-law mass distribution P(m) ~ jT^-TconBcrvcd^ 
This distribution coexists with an infinite aggregate which accommodates a finite fraction of the total mass — a real 
space analog of Bose-Einstein condensation p3| . The exponent T conscrvc d was found to be 5/2 within mean field theory 
p3| , p4[ and ~ 2.33 in 1-d p3| , and the same exponent was found to describe P(m) at the critical point. Evidently, 
the lack of mass conservation in the In-out model is responsible for the absence of the infinite aggregate in its high q 
phase, as well as the change in the power to r in the Takayasu phase and r c at criticality. 
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Another interesting difference between the conserved and the In-out modei is the effect of a preferred direction for 
the motion of masses (a mass at site i hops with a higher probability to i — 1 than i + 1). We have checked that such 
a bias does not change the critical exponents of the In-out model. However, for the conserved mass model, the bias 
in direction changes the value of the exponents at the transition and in the aggregate phase [ p3[ . 

The phase transition in the In-out model has some interesting implications for noncquilibrium wetting transitions if 
we interpret the configuration of masses as an interface profile regarding m.j as a local height variable. Although the 
dynamics of the mass profile in our In-out model is not physical when interpreted as interface dynamics, nevertheless 
the phase transition in our model can be qualitatively interpreted as a nonequilibrium wetting transition of the 
interface. In the In-out model, the fact that the mass at each site is necessarily non-negative translates into the 
restriction that the mass profile is always above a wall at a fixed height (in our case 0). The presence of this 
constraint is the key factor for the wetting transition. At fixed p, as we increase q, the mean height (to) does not 
grow with time as long as q < q c . This is our 'Exponential' phase where the mass or the interface profile is bound 
to the substrate at zero height. This phase is also 'smooth' as the mean square height fluctuation does not grow 
with system size. For q > q c , the interface unbinds from the substrate and the mean height (to) ~ vt grows linearly 
with time with a velocity v. This phenomenon is similar to 'wetting' or 'depinning' of interfaces in general. In this 
'wet' phase (Takayasu phase in our model), the interface is rough. Unlike recently studied models of nonequilibrium 
wetting, where the interface in the growing phase is self-affine pl| , ^9| , our model describes a much rougher interface 
for q > q c . At the transition q = q c , though, the interface is self-affine with a roughness exponent \ — Z P — 0-7- 

There are various open questions that remain to be settled. In this paper, we have only studied the phase transition 
in the steady-state single site mass distribution function. It would be very interesting to study the spatial correlations 
between masses at different sites and to track the behaviour of mass-mass correlation function as one crosses the phase 
boundary in the p — q plane. 

Also in this paper we have only studied the simplest model where the rates of adsorption, desorption and hopping 
are constants and independent of particle mass. An important question is whether this phase transition would persist 
for general mass-dependent rates. In earlier work a model with aggregation, adsorption and desorption was 
studied, but no transition to a power-law phase was found; the difference is traceable to the fact that in that model, 
the rate of removal of mass is proportional to the mass, unlike the unit-mass desorption process considered in the 
In-out model. It is therefore highly desirable to identify the class of models with mass-dependent rates where the 
phase transition described here will persist. 
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a? tor 
q 

FIG. 1. The function s(q) as a function of q for p = 1 (shown by the solid line) within the mean field theory. It deviates 
from the dotted line (s(q) = g) for q > q c = 3 — 2\/2- 
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FIG. 2. The phase diagram of the In-out model in the p — q plane. The dotted line denotes the mean field phase boundary 
q c (p) = p + 2 — 2^/p + 1 and the triangles mark the numerically obtained critical points in 1-d. 
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FIG. 3. The steady-state mass distribution P(m) vs. m for the In-out model in 1-d. The value of q is kept fixed at q = 1 
and the data is shown for three representative valuesof p, respectively p <, =, > p c ~ 2.35. 




FIG. 4. The velocity v as a function of q for fixed p = 2.35 in 1-d. The velocity is zero for q < q c 
(q — q c ) y for q > q c with y ~ 1.47 in d = 1. 



1.0 and increases as 



10 




Q 2 I I I I I I I I 

'0 1 2 3 4 5 6 7 8 

-2 

tL 

FIG. 5. Scaling plots of the finite-size studies of width versus time for four different system sizes L = 16, L = 32, L = 64 
and L = 128 at the critical point p c « 2.352 with fixed q = 1. The width is expected to follow the scaling form W w t 13 Z(t/L z ). 
The best collapse is obtained for (3 ~ 0.358 with 2 = 2. 
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